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Abstract: The focus of our research was to automatically generate the most adequate vegetation density maps for 
orienteering purpose. Application Karttapullatuin was used for automated generation of vegetation density maps, 
which requires LiDAR data to process an automatically generated map. A part of the orienteering map in the area of 
Kazlje-Tomaj was used to compare the graphical display of vegetation density. With different settings of parameters 
in the Karttapullautin application we changed the way how vegetation density of automatically generated map was 
presented, and tried to match it as much as possible with the orienteering map of Kazlje-Tomaj. Comparing more 
created maps of vegetation density the most suitable parameter settings to automatically generate maps on other are- 
as were proposed, too. 
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1. Introduction 

Orienteering maps, maps for orienteering sport, are the most standardized type of the maps (Zentai, 2000), de- 
fined by International Specification for orienteering Maps (ISOM). Standardization is essential to assure equal and 
fair conditions for all competitors. It is also very important, that the maps are regularly updated, showing the actual 
situation on the field. Maps are very detailed and present landscape at large scale, which allows competitor to esti- 
mate the time needed to move along the course. Besides the other criteria influencing speed, like steepness, ground 
roughness, physical built obstacles, etc. vegetation is probably the one causing the highest running speed reduction. 
That’s also the reason, why vegetation is one of the most important elements of orienteering map; there are even 22 
of all 121 symbols and two different colors of all together five used for vegetation features: green, yellow and addi- 
tionally also white, uncolored areas, which also shows the clear, fast runnable forest area. Figure 1 shows different 
vegetation areas and objects on two inserts of orienteering map. 



Fig. 1 . Different symbols for vegetation on Orienteering map (maps by Ivan Nagy) 
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On the other hand, there are many problems of presenting vegetation on orienteering maps: lack of useful data- 
bases and other sources of data (comparing to other content, like buildings, roads), vegetation type and density 
changes in small areas and in short period and finally, interpretation of the density, height and type of vegetation and 
their influence on run speed reduction. Such interpretation in most cases can be done only by the experienced map 
maker, former orienteering runner in most cases; on the terrain only. 


2. LiDAR data as the source for vegetation density map 


LiDAR technology was recognized as very useful as the main source data for creating orienteering and also topo- 
graphic maps (Ditz et al, 2014, Petrovic 2010, 201 1, Petrovic and Podobnikar 2011). The LiDAR point cloud as the 
basic result of laser scanning includes reflections of different terrain objects, from ground, built objects, and also 
grass, leaves, branches and other parts of vegetation (Figure 2). The resolution of cloud points can be very high and 
therefore representation of the terrain with all objects on it may be quite detailed. 
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Fig. 2. Point cloud shows different types of vegetation, grassland, bushes, and trees 

Such point cloud is probably the best possible source data for orienteering map maker to get the feeling about the 
density and sometimes even the type of vegetation, even better than camera movie, which shows only two dimen- 
sional, flat images. From point cloud map maker have to estimate the run speed reduction due to the vegetation 
height, density and type. Different sets of classification, using intensity, order of return, and height above ground 
etc. can support him to understand and recognize complex vegetation appearance. But even from those classified da- 
ta to vegetation content on the orienteering map there is still intensive manual work. 

Some applications try to automatize the process of interpretation for different content from lidar point clouds, 
where exact purpose of derived data has to be predefined. Orienteering maps with global standardizations are quite 
appropriate for that. One such freely available application Karttapullautin can be used for automatic creation of ori- 
enteering maps with all content. Maps are not perfect, but good enough for trainings or small local events (Petrovic 
2013, 2014). One of the strength of the application is the parameter tuning possibility; every map creator can adjust 
number of parameters according to his specifics: density and complexity of lidar point cloud data, type of the terrain, 
including type and density of vegetation, type of ground (rocky, sandy, moody) etc. More than 30 parameters can be 
adjusted for vegetation presentation only and all of them partly influence the final result. 
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3. Automatically created vegetation areas for orienteering map 


The aim of our research was therefore finding the optimal set of parameters for selected specific area (Cerar, 
2016) using application Karttapullautin. As a reference we used the existing orienteering map, made by experienced 
map maker, who mapped the vegetation based on field check. The whole methodology based on iterative tuning sets 
of parameters with comparing results to reference after each step. The success of result was calculated with the over- 
lapping ratio of the same and neighboring category between reference and the result. ISOM defines 13 different are- 
al types of vegetation: from open (grass) land to impassable dense vegetation area, from pure natural to planted cul- 
tivated area. In our research we focused on five most important for the runner: 

(1) open land (full yellow), 

(2) open runnable forest (white), 

(3) forest that enables slow running (light green), 

(4) forest where running is difficult (middle green) and 

(5) impassable forest or other vegetation (dark green). 

Automatically created vegetation map was compared with the reference map, we measured the areas on both 
maps covered with the same category and marked them as a full coverage (F). But usability of the presented vegeta- 
tion can’t be estimated only with an exact overlapping. Sometimes it is very difficult to differ between neighbor cat- 
egories on the terrain, but middle green would give the similar information to the runner as dark green, although not 
the same. Therefore we decided to count also the areas overlapping with neighboring category (N), but weighted. 
The final usability (U) of the vegetation on the map was therefore calculated as: 

U = (F + 0,65 x N) / 10, 

where F is % of exact overlaping of the category in the selected area (Full coverage) and 

N is % of overlapping neighboring category, as presented in Figure 3 



(1) open land (full yellow), 

(2) open runnable forest (white), 

(3) forest enables slow running (light green), 

(4) forest where running is difficult (middle green) and 

(5) impassable forest or other vegetation (dark green). 


F exact overlap of the category 


N overlap neighboring category 


<3 automatically created (3) light 


green 


Fig. 3. Calculation of the usability of vegetation map 
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3.1 Case area 

Area in Karst region in Slovenia was used as a case area. Karst area is very popular for orienteering in Slovenia, 
karstic plateau with complex relief formations (depressions, pits, ditches), stony ground and objects (boulders, rock 
faces, stone walls) and diverse vegetation attracts athletes from many European countries for trainings and competi- 
tions. Due to mild climate running is possible from early spring to late autumn, usually also in winter, when lot of 
other terrains in central Europe are not suitable for orienteering. Therefore a lot of maps that are regular updated ex- 
ist on the area and one of those maps, Kazlje - Tomaj (Fig. 4) was used as a reference in our research. 



Fig. 4. Case area in Karst region in Slovenia, reference map Kazlje - Tomaj (map by Ivan Nagy) 


3.2 Used data 

Slovenia is one of the rare countries, where whole territory is covered with high density laser scanning data and 
where all the data is freely available. The project of LiDAR scanning of Slovenia started in 201 1 as a part of flood 
protection measures. Partly in 201 1, but mostly in 2014 and 2015 whole territory was surveyed with density 5 points 
/ sq. m. Web service at National Agency for the Environment offer users data in different formats (classified point 
cloud, point cloud of ground points and Dim DTM) in two different coordinate systems (Fig 5). 


udar“ 



Fig. 5. Web service with LiDAR data for Slovenia (http://gis.arso.gov.si/evode/) 
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Data is available in 1 sq. km blocks, for our research we needed complete point cloud with ground and other 
points, too. 


3.3 Case study 


Kartapullautin is free application for noncommercial use or if used for navigation sport mapping (orienteering, 
rogaining, adventure racing map making), created by Finnish author Jarkko Ryyppo 
( http://www.routegadget.net/karttapullautin/ ). It can be used for generating orienteering training maps without field 
checking. In our research we used version 20150426; the latest version 20160405 was issued just in the time of our 
research. Point cloud LAS data is used as a source data while the final result is image of complete orienteering map, 
composed of partly created vector data (contours, point objects) and set of raster images (buildings, water areas, 
vegetation). Vegetation mask as one separate raster layer can therefore be analyzed separately from other content. 39 
parameters that can be adjusted by user in Kartapullautin and influence on final vegetation mark image are intro- 
duced in Table 1 with initial predefined values. 


Table 1 . Parameters for vegetation mask in Kartapullautin 


greenground 

0.8 

zonel 

0.5 1 4 1 99 1 1 

zone2 

4 1 7 1 10.0 1 0.4 

thresholdl 

0.20| 2.0 10.45 

threshold2 

2.0 1 6 10.41 

threshold3 

6 1 16 10.42 

threshold4 

16 1 20 10.46 

thresholds 

20 1 99 10.47 

pointvolume factor 

0.35 

firstandlastreturnfactor 

0.1 

lastreturnfactor 

0.2 

greenshades 

0.4 1 0.8 1 1.3 1 2.6 1 99 

lightgreentone 

230 

greendetectsize 

3 

groundboxsize 

3 

medianboxsize 

12 

medianboxsize2 

1 

yellowheight 

0,8 

yellowthresold 

0,89 


Greenshades define range of different used green colors, from light to dark green. Threshold 1-5 parameters are 
used for balancing tall and low forest. Very low vegetation mostly influence the run speed, middle one visibility, 
while tall vegetation only define type of the forest. Zone 1 and 2 parameters are used for separating different types 
of forest trees, according to density of lover branches. The other parameters are used to divide areas between open, 
semi open, rough and forest areas and to determine the smallest area that is presented as a separate area. Too small 
recognized and presented areas would lead to extremely scattered and therefore poorly readable map. 
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Using initially set parameters we created the first version of automatically generated vegetation mask map and 
compered it with the reference map (Fig. 6). 



500 m 


Fig. 6. Automatically generated vegetation mask; reference map; overlapping analyze (pink shows full coverages areas and blue over- 
lapping neighboring category) (Cerar, 2016) 

Calculation according to presented methodology gave results presented in Table 2. 


Table 2. Results of analyze using initial parameters 



area on autom. 
gener. map (sq. m) 

(F) % of ex- 
act overlap 

(N) % of overlapping 
neighbouring category 

(U) usa- 
bility 

(1) open land (full yellow), 

89634,52 

83,62 

7,60 

8,86 

(2) open runnable forest (white), 

69953,79 

51,16 

39,54 

7,69 

(3) forest enables slow running (light green), 

181995,55 

28,02 

45,58 

5,76 

(4) forest where running is difficult (middle green) and 

282444,97 

22,79 

43,27 

5,09 

(5) impassable forest or other vegetation (dark green). 

90375,01 

48,51 

27,54 

6,64 

not classified 

11276,16 




sum 

725680,00 

38,02 

37,26 

6,22 

average 

142880,77 

46,82 

32,71 

6,81 


Initial values of parameters gave partly satisfied results but our goal was to see, how far usability of the vegeta- 
tion mask map can be improved with adjusting parameters. The iterative procedure was used, starting with the first 
group of parameters, continuing with next ones. Figure 7 shows influence of adjusting parameter “greenshades” on 
final vegetation mask map. 
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Fig. 7. Influence of adjusting parameter “greenshades” on vegetation mask map (Cerar, 2016) 

After 4 complete iterations where in each most of parameters were adjusted we managed to improve the overall 
usability in 9 % for entire area and in 2 % only in average of all vegetation categories, as presented in Table 3. Fig- 
ure 8 shows the overlapping map after fourth iteration and Table 4 the final values for parameters. 
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Table 3. Results of analyze after fourth iteration 



area on autom. 
gener. map (sq. m) 

(F) % of ex- 
act overlap 

(N) % of overlapping 
neighbouring category 

(U) usa- 
bility 

(1) open land (full yellow), 

90224,79 

83,08 

7,48 

8,79 

(2) open runnable forest (white), 

203059,09 

40,83 

45,79 

7,06 

(3) forest enables slow running (light green), 

209971,89 

30,88 

44,88 

6,01 

(4) forest where running is difficult (middle green) and 

120367,90 

30,02 

46,63 

6,03 

(5) impassable forest or other vegetation (dark green). 

92056,32 

50,80 

27,67 

6,88 

not classified 

10000,01 




sum 

725680,00 

42,70 

38,50 

6,77 

average 

143136,00 

47,12 

34,49 

6,95 


Table 4. Final values of parameters for vegetation mask in Kartapullautin 


greenground 

0.8 

zonel 

0.5 14 1 99 1 1 

zone2 

4| 7 1 10.0 1 0.4 

thresholdl 

0.2012.5/0.28 

threshold2 

2.5/5/0.37 

thresholds 

5/8/0.43 

threshold4 

8/12/0.39 

thresholds 

12/99/0.37 

pointvolumefactor 

0.15 

firstandlastreturnfactor 

0.1 

lastreturnfactor 

0.2 

greenshades 

0.64|0.98|1.37 

lightgreentone 

205 

greendetectsize 

3 

groundhoxsize 

3 

medianboxsize 

12 

medianboxsize2 

1 

yellowheight 

0,88 

yellowthresold 

0,89 
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Fig. 8. Overlapping analyze with initial and final values of parameters (pink shows full coverages areas and blue overlapping neighbor- 
ing category) (Cerar, 2016) 


3.4 Testing parameters in other areas 

Although adjustments of parameters only slightly improved the results compared with initial parameters we tried 
to verify, what happens in other areas. Since all adjustments were done for karstic area we select another area in 
Karts plateau, few km from initial case area. Results are presented in Fig. 9 and Table 5. 



Fig. 9. Reference orienteering map for area 2 and overlapping analyze with final values of parameters (pink shows full coverages areas 
and blue overlapping neighboring category) (Cerar, 2016) 


Table 5. Usability of vegetation mask map in area 1 and 2 



(U) usability in area 1 

(U) usability in area 2 

(1) open land (full yellow), 

8,79 

9,39 

(2) open runnable forest (white), 

7,06 

7,03 

(3) forest enables slow running (light green), 

6 

5,8 

(4) forest where running is difficult (middle green) and 

6,03 

2,98 

(5) impassable forest or other vegetation (dark green). 

6,88 

3,96 

whole map 

6,77 

7,19 
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4. Conclusions 

The final test in area 2 showed even better results of overlapping vegetation areas as in the area, where parame- 
ters were adjusted. This might point out that adjustment in the first test area wasn’t optimal, while different vegeta- 
tion combination might also be a reason for such result. Further on, another two tests with same parameters were 
made in two completely different areas with different vegetation, where results were significantly worse, as ex- 
pected. More tests using different methods and different applications should be done before finding optimal solution 
for creating vegetation map mask for orienteering maps. 
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